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Abstract - We study the sample paths properties of Operator scaling Gaus- 
sian random fields. Such fields are anisotropic generalizations of anisotropic 
self-similar random fields as anisotropic Fractional Brownian Motion. Some 
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these anisotropies and thus to the geometry of these fields. 
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1 Introduction and motivations 

Random fields are now used for modeling in a wide range of scientific areas 
including physics, engineering, hydrology, biology, economics and finance 
(see [30] and its bibliography). An important requirement is that the data 
thus modelled present strong anisotropies which therefore have to be present 
in the model. Many anisotropic random fields have therefore been proposed 
as natural models in various areas such as image processing, hydrology, 
geostatistics and spatial statistics (see, for example, Davies and Hall [15], 
Bonami and Estrade |8], Benson et al. [3]). Let us also quote the exam- 
ple of Levy random fields, deeply studied by Durand and Jaffard (see [H]), 
which is the only known model of anisotropic multifractal random field. In 
many cases, Gaussian models have turned to be relevant when investigating 
anisotropic problems. For example the stochastic model of surface waves 
is usually assumed to be Gaussian and is surprisingly accurate (see [2T]). 
More generally anisotropic Gaussian random fields are involved in many 
others concrete situations and then arise naturally in stochastic partial dif- 
ferential equations (see, e.g., Dalang [M], Mueller and Tribe [23], Oksendal 
and Zhang [26], Nualart [25j). 

In many situations, the data present invariant features across the scales 
(see for example ^). These two requirements (anisotropy and self-similarity) 
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may seem contradictory, since the classical notion of self-similarity defined 
for a random field {X{x)}^^^d on by 

{X{ax)}^^^d = {a^oX(x)}^giRd , (1.1) 

for some Hq G M (called the Hurst index) is by construction isotropic and 
has then to be changed in order to fit anisotropic situations. To this end, 
several extensions of self-similarity property in an anisotropic setting have 
been proposed. In [19] , Hudson and Mason defined operator self-similar pro- 
cesses {X(i)}fgK with values in W^. In |20j, Kamont introduced Fractional 
Brownian Sheets which satisfies different scaling properties according to the 
coordinate axes. More recently, in [7] Bierme, Meerschaert and Scheffler 
introduced the notion of Operator Scaling Random Fields (OSRF). These 
fields satisfy the following anisotropic scaling relation : 

{X(a^»:E)},gK. ^ {a^»X(x)},gK<^ , (1.2) 

for some matrix Eq (called an exponent or an anisotropy of the field) whose 
eigenvalues have a positive real part and some Hq > (called an Hurst index 
of the field). The usual notion of self-similarity is extended replacing usual 
scaling, (corresponding to the case = Id) by a linear scaling involving 
the matrix Eq (see figure [1] below). It allows to define new classes of random 
fields with new geometry and structure. 




Ag{i,---,io} Ag{i,---,io} Ae{i,---,io} 

Figure 1: 

Action of a linear scaling x i— )■ \^x on the smallest ellipsis. 

This new class of random fields have been introduced in order to model 
various phenomena such as fracture surfaces (see [27] ) or sedimentary aquifers 
(see [1]). In [7], the authors construct a large class of Operator Scaling Sta- 
ble Random Fields with stationary increments presenting both a moving 
average and an harmonizable representation of these fields. 
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In order to use such models in practice, the first problem is to recover 
the parameters and Eq from the inspection of one sample paths. Even if 
we consider the model mentioned above in the Gaussian case, the problem 
of identification of an exponent of self-similarity Eq (which in some case is 
not unique) and of an Hurst index Hq is an open problem. 

The first step in the resolution of this question involves an identification 
of some specific features of exponents and indices which can be recovered 
on sample paths. This paper is a first step : we will prove that from the 
regularity point of view these exponents and Hurst indices satisfy what we 
call optimality properties. More precisely, we prove that (see Theorem 14. ip . 
the Hurst index Hq maximizes the local critical exponent of the field in 
specific functional spaces related with the anisotropy matrix Eq among all 
possible critical exponents in general anisotropic functional spaces. 

Therefore, the results of the present paper open the way to the following 
strategy to recover the Hurst index. One first have to consider a discretized 
version of the set of all possible anisotropics. In each case an estimator of 
the critical exponent related with these anisotropics has to be given. There- 
fore, one has to locate the maximum of all these estimators-which can be 
based on anisotropic quadratic variations-and to identify the corresponding 
values of the anisotropy. Combining the resluts of this paper and these of 
m, our approach can thus be turned in an effective algorithm for the esti- 
mation of the anisotropy of self similar textures (see [281 129]). The study of 
the estimators related to critical indices in anisotropic Besov spaces from a 
statistical point of view will be the purpose of a forthcoming paper. 

Our optimality result comes from sample paths properties of the model 
under study in an anisotropic setting. This approach is natural : In [20| , Ka- 
mont studied the regularity of the sample paths of the well-known anisotropic 
Fractional Brownian Sheet in anisotropic Holder spaces related to Fractional 
Brownian Sheet. Moreover, some results of regularity in specific anisotropic 
Holder spaces related to matrix Eq have already be established for operator 
scaling self-similar random fields (which may be not Gaussian) in [6] or in 
the more general setting of strongly non deterministic anisotropic Gaussian 
fields in [37] . We then extend already existing results by measuring smooth- 
ness in general anisotropic spaces not necessarily related to the exponent 
matrix Eq of the field. 

This paper is organized as follows. In Section [21 we briefiy recall some 
facts about Operator Scaling Random Gaussian Fields (OSRGF) and de- 
scribe the construction of [7J of the model. In Section [3l we present the 
different concepts used for measuring smoothness in an anisotropic setting 
and especially anisotropic Besov spaces. Section |4] is devoted to the state- 
ment of our optimality and regularity results. Finally, Section [5] contains 
proofs of the results stated in Section jH 
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In the sequel, we will use some notations. For any matrix M 

Pmm(Af) = min (|Re(A)|), p^UM) = ^ max (|Re(A)|) , 

AeSp(M) xeSp{M) 

where Sp(M) denotes the spectrum of matrix M. 
For any real a > 0, denotes the matrix 



a^' = exp(Mlog(a)) = ^ ^^-^ 



A:>0 



In the following pages, we denote f the collection of matrices of 
whose eigenvalues have positive real part. 



2 Presentation of the studied model 

The existence of operator scaling stable random fields, that is random fields 
satisfying relationship (jl.2p . is proved in [7]. The following theorem (Theo- 
rem 4.1 and Corollary 4.2 of [7J) completes this result by yielding a practical 
way to construct a Operator Scaling Stable Random Field (OSRF) with sta- 
tionary increments for any Eq £ and Hq G (0, Pmin(-E'o))- We state it 
only in the Gaussian case, having in mind the problem of the estimation of 
the Hurst index Hq and the anisotropy Eq. 

Theorem 2.1 Let Eq he in £^ and p a continuous function with positive 
values such that for all x ^ 0, p{x) ^ 0. Assume that p is ^EQ-homogeneous, 
that is : 

Va > 0, G M'^, p(a^o^) = ap{^) . 
Then the Gaussian field 

X,ix) = [ (e'<-'«> - l)p(0-^"-^dW(0 , (2.3) 

exists and is stochastically continuous if and only if Hq G (0, Pmin(-£'o))- 
Moreover this field has the following properties : 

1. Stationary increments : 

yh G {X,ix + h)- X,{h)},^^, ^4'^ {X,{x)},^^, . 

2. The operator-scaling relation \1.S\) is satisfied. 

Remark 2.1 The assumption of homogeneity on the function p is necessary 
to recover linear self-similarity properties of the Gaussian field {Xp{x)}^^^d. 
The assumption of continuity on p ensures that the constructed field is 
stochastically continuous. 
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Remark 2.2 In general, the couple {Ho,Eq) of an OSRF is not unique. 
Indeed, if Hq and Eq are respectively an Hurst index and an exponent of 
the OSRF {X{x)}^fz^d, then for any A > so do XHq and XEq. 
Uniqueness of the Hurst index Hq can be recovered by choosing a normaliza- 
tion for £^01 for example Tr(£'o) = d. However, even under this assumption, 
Eq is not necessarily unique. Nevertheless remark that, under the assump- 
tion Tr(£'o) = d, two anisotropies of an OSRF have necessarily the same real 
diagonalizable part (see Section 15.21 for a definition). We refer to Remark 
2.10 of [7J for more details on the structure of the set of exponents of an 
OSRF. 

Remark that Theorem 12.11 relies on the existence of ^E'Q-homogeneous 
functions. Constructions of such functions have been proposed in [7] via an 
integral formula (Theorem 2.11). An alternative construction, more fitted 
for numerical simulations, can be found in jl3j . 

3 Anisotropic concepts of smoothness 

Our main goal here is to study the sample paths properties of this class of 
Gaussian fields in suitable anisotropic functional spaces. This approach is 
quite natural (see [201 E]) since the studied model is anisotropic. To this 
end, suitable concepts of anisotropic smoothness are needed. The aim of 
this section is to give some background about the appropriate anisotropic 
functional spaces : Anisotropic Besov spaces. These spaces generalize clas- 
sical (isotropic) Besov spaces and have been studied in parallel with them 
(see [U [lOj for a complete account on the results presented in this section) . 
The definition of anisotropic Besov spaces is based on the concept of pseudo- 
norm. We first recall some well-known facts about pseudo-norms which can 
be found with more details in |22j . 

3.1 Preliminary results about pseudo-norms 

In order to introduce anisotropic functional spaces, an anisotropic topology 
on is needed. We need to introduce a slight variant of the notion of 
pseudo-norm introduced in [22], fitted to the case of discrete dilatations. 

Definition 3.1 Let E G A function p defined on is a (M'^,£') 

pseudo-norm if it satisfies the three following properties : 

1. p is continuous on W^, 

2. p is E -homogeneous, i.e. p{a^x) = ap{x) Vx € W^, Va > 0, 

3. p is strictly positive on \ {0} . 
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For any (R'^,£') pseudo-norm, define the anisotropic sphere Sq{p) as 
S§{p) = {xe M^; p{x) = 1} . (3.4) 




Figure 2: 

Examples of anisotropic spheres for different anisotropies. 

Proposition 3.1 For all x € M'^\{0}, there exists an unique couple {r,9) G 

^ ^oiP) -^^c/i that X = r^B. 

Moreover Sq{p) is a compact o/M'* and the map 

(r, e)^x = r^e , 

is an homeomorphism from M.\ x Sq{p) to \ {0}. 

The term "pseudo-norm" is justified by tlie following proposition : 

Proposition 3.2 Let p a (W^jE) pseudo-norm. There exists a constant 
C > such that 

p{x + y)<C{p{x) + p{y)), Vx, yGM'^. (3.5) 

The following key property allows to define an anisotropic topology on R'^ 
based on pseudo-norms and then anisotropic functional spaces : 

Proposition 3.3 Let p\ and p2 he two {W^, E) pseudo-norms. They are 
equivalent in the following sense : There exists a constant C > such that 

^Pi{x) < P2{x) < Cpi{x), yxeM.'^ . 

In particular, two different {W^, E) pseudo-norms define the same topology 
on W^. 
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3.2 Anisotropic Besov spaces 

Let E e and | • [tg a fixed (M'^,*E')-pseudo-norm. For xq G M*^ and r > 0, 
B|.|^^(a;o, r) denotes the anisotropic ball of center xq and radius r, namely 

5|.|^^(xo,r) = {x€ R'^,\x - xo\tE< r} . 



Definition 3.2 Let i/j^ G <S(M^) 6e swc/i i/tai 

I v?(e) = o*/ici*E>2. 

For any positive integer j, set 

vf(o = vl(2-^''''0-v?(2-(^-n). 

Then 

is an anisotropic partition of the unity satisfying supp(^J') C S|.|^^(0, 2-''^-'-)\ 

The anisotropic Besov spaces Bp q{'R'^, E) are then defined as follows : 
Definition 3.3 Let < p,q < oo and s G M. Define 

II/IIb|.M= . (3.6) 

Then 

b;^^{r'',e) = {/ g ^'(m'^), WfWB^^^iR^E) < +^} ■ 

The matrix E is called the anisotropy of the Besov space Bp^g{W^, E). 
In a more general way, i/ /? G M, define 

\j=o 

Then 
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Remark 3.1 One can prove that this definition is independent of the choice 
of the function ipQ involved in the definition of the Besov space q{W^, E). 

Remark 3.2 Let E ^ £^ and | • 1 1^; a (M'^, pseudo-norm. For any A > 0, 
I • is a (M'^,A*^) pseudo-norm. Hence for any s > 0, B^^g{R'^,XE) = 
B;,g{^'',E). 

So, without loss of generality, we assume in the sequel that Tr(£') = d. 
We then define 

4 = {Ee Tt{E) = d} , 

where is the collection of d x d matrices whose eigenvalues have positive 
real part. As it is the case for isotropic spaces, anisotropic Holder spaces 
C'^{M.'^,E) can be defined as particular anisotropic Besov spaces. 

Definition 3.4 Let s be in M and /? € M. The anisotropic Holder spaces 
C'(R'^,E) andC^^^^^,^{M.'^,E) are defined by 

C^{R'',E) = <^rid q,^^^,{R^,E) = B'^^^^^,^^^,{R^, E) . 

Proposition 3.4 Let < s < PminiE) and /3 G M. Then the two norms 

, \f{x + h)-f{x)\ 

Loo(Kd) + sup sup 



h\E<ixmd V l/illllogd^b)!^ 



defined by / I A'. 6]) . 

oo,oo,| log \ ^ 

are equivalent in C^^^^^p{R'^, E). 



Remark 3.3 Anisotropic Holder spaces admit a characterization by finite 
differences of order M > 1 under the general assumption s > 0. Here, 
we only need to deal with the case < s < pram{E) and have thus stated 
Proposition 13.41 in this special setting. 



Let us comment Proposition 13. 4[ Let < s < Pmin 

(E) and € M. A 

bounded function / belongs to C^^^^^p{M.'^, E) if and only if for any r G (0, 1), 
e 5^(1 ■ l^;) and x G W^, one has 

|/(x + r^G)-/(x)| <Cor^|log(r)|^ , 

for some Cq > 0. 

Hence, a function / belongs to the Holder space q, |;3(M'^,£;) if and only if 
its restriction /e along any parametric curve of the form 

r > r^G , 
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with e G S^{\ ■ \e) is in the usual Holder space Cf^ and H/ellc (r) 

does not depend on B. Roughly speaking, the anisotropic "directional" 
regularity in any anisotropic "direction" has to be larger than s. In other 
words, we replace straight lines of isotropic setting by curves with parametric 
equation r > i-7> r^Q adapted to anisotropic setting. 

To state our optimality results we need a local version of anisotropic Besov 
spaces : 

Definition 3.5 Let E € £^ be a fixed anisotropy, < p,q < oo, < s < oo 

and f GLfjR'i). 

The function f belongs to Bp^i^^{W^,E) if for any € P(M'^), the function 
iff belongs to B^^g{R'^,E). 

The spaces B"^^ \ \og\f loc^^'^^^^ defined in an analogous way for any 

< p,q < oo, < s < oo, (3 eR. 

The anisotropic local critical exponent in anisotropic Besov spaces i?* ^(M'^, 
of f e L^^^(IR'^) is then defined by 

afjoc(.E,p,q) = sup{s, / G ,^ ;„^(M'^, S)} . 

In the special case p = q = oo, this exponent is also called the anisotropic 
local critical exponent in anisotropic Holder spaces of / € L^^(M'^) and is 
denoted by a/,ioc 

4 Statement of our results 

In what follows, we are given Eq € and a (M'^, *Eq ) pseudo-norm. 
We denote {Xp^ ,Hoix)}xGR'^ OSRGF with exponent Eq and Hurst index 
Hq defined by (|2.'3p with p = pE^- 

We first state our optimality result and characterize in some sense an 
anisotropy Eq and an Hurst index of the field {Xp^^^Ho{x)}x^M.d. These re- 
sults come from an accurate study of sample paths properties of the OSRGF 
{Xp^^^Hoix)}x^^d in anisotropic Besov spaces (see in Section ??). 

We assume -without loss of generality -that Eq G £'j", namely that 
all the eigenvalues of Eq have a positive real part and that Tr(So) = d. 
Our results will be based on a comparison between the topology related 
to the pseudo-norm peq involved in the construction of the Gaussian field 
{Xp^^^Ho{x)}xeR'i and this of the analyzing spaces Bp^g{W^, E). To be able 
to compare these two topologies, we also assume that E ^ £^ . 

The main result of this paper is the following one : 
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Theorem 4.1 Let {p,q) G [l,+oo]'^ and Eq E SJ^ . Then almost surely 
axpj^ jj^,loc{Eo,p,q) = sup{aXp^ jj^^ioc{E,p, q), E G £j,E commuting with Eq} 



=-0 

= Ho , 

that is the value E = Eo maximizes the anisotropic local critical exponent of 
the OSRGF {Xp^^^Hoix)}x^K'i among all possible anisotropic local critical 
exponent in anisotropic Besov spaces with an anisotropy E commuting with 

Eq. 

Remark 4.1 Since E and Eo are commuting, these matrices admit the 
same spectral decomposition. Hence, in fact we proved that any anisotropy 
£^0 maximize the critical exponent among matrices having the same spectral 
decomposition. Thus, in the general case, we implicitly assumed that the 
spectral decomposition of anisotropy matrix is known. In dimension two, 
we have a stronger optimality result about anisotropy £"0 and Hurst index 
Ho, involving matrices of £^ which do not commute necessilary. 

To prove Theorem 14. 1^ we investigate the local regularity of the sample 
paths of {Xp^^^^Ho{x)}^^-^d in general anisotropic Besov spaces. But before 
any statement, we first need some background about the concept of real 
diagonalizable part of a square matrix. This notion is based on real additive 
Jordan decomposition of a square matrix (see for e.g. to Lemma 7.1 chap 9 
of [18] where a multiplicative version of Proposition 14.11 is given) : 

Proposition 4.1 Any matrix M of M^iM.) can he decomposed into a sum 
of three commuting real matrices 

M = D + S + N , 

where D is a diagonalizable matrix in Mrf(M), S is a diagonalizable matrix 
in Mrf(C) with zero or imaginary complex eigenvalues, and N is a nilpotent 
matrix. Matrix D is called the real diagonalizable part of M, S its imaginary 
semi-simple part, and N its nilpotent part. 

Now we are given two commuting matrices Eo, E of Sj^ . Let Do (resp 
D) be the real diagonalizable part of matrix Eq (resp E). Since matrices 
Eq and E are commuting, so do matrices Do and D. Furthermore, matrices 
Do and D are diagonalizable in M(i{^) then they are simultaneously diago- 
nalizable. Up to a change of basis, we may assume that Dq and D are two 
diagonal matrices. More precisely, suppose that 



D, 







V 




D 



V 




(4.7) 
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with 

^<-<^. (4.8) 

Since Tr(ii^o) = Ti{E) = d, one has Xm/^m ^ 1- 

The regularity results about sample path of the field {Xp^^^Ho{x)}xeR'' 
are summed up in the following theorem. 

Theorem 4.2 Let 1 < p < +oo, I < q < +oo. Almost surely the anisotropic 
local critical exponent dXp^ u^Mci^^P^ l) anisotropic Besov spaces ^(M'^, E) 
of the OSRGF {Xp^^^Ho{x)}^^^d satisfies 

ax HaME,p,q) = — ^ < Hq . 

In particular, in the special case E = Eq, one has ctXp^ HQ,iociE,p, q) = Hq. 

In other words Theorem 14.21 asserts that when one measures local regularity 
of the sample paths along anisotropic directions different from those asso- 
ciated to an anisotropy of the field Eq, one loses smoothness. The further 
the anisotropic direction of measure from the genuine anisotropic direction 
associated to the field are, the smaller the anisotropic local critical exponent 
is. This anisotropic local critical exponent can take any value in the range 
iO,Ho]. 

The special case p = q = +oo yields us the following result about 
anisotropic Holderian regularity of the sample paths. 

Corollary 4.1 Almost surely the anisotropic local critical exponent of the 
sample paths o/{Xp^^^^j^(,(x)}^g]gd in anisotropic Holder spaces equals {XmHQ)/X^ 
and is always lower than Hq . In particular, if E = Eq this critical exponent 
equals the Hurst index Hq. 

Remark 4.2 This estimate on anisotropic local critical exponent was al- 
ready known in the case E = Eq (see [6]). 

Theorem 14.21 allows us to obtain regularity results which extend those 
proved in the case j? = g = oo in the usual isotropic setting. Since matrices 
£"0 and Id are commuting, we can apply the above result to the case E = Id. 
Note that in this case = /Oniax(-E'o)- We obtain the following proposition 



Proposition 4.2 Almost surely the local critical exponent of the sample 
paths of {Xp^^^^Ho{'^)}xeM.'i classical Besov spaces equals -ffo/Pmax(-So). 
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In particular, for p = q = oo, almost surely the local critical exponent 
of the sample paths of {Xp^^^Ho{x)}x£R<i 'in classical Holder spaces equals 

Ho/pmaxi^o)- 

Remark 4.3 In the special case p = q = oo, we recover already known 
results about classic Holderian regularity (see Theorem 5.4 of [7]). Re- 
call that this theorem is based on directional regularity results about the 
Gaussian field {Xp^^^^Ho} comes from an estimate of the variogram 
vx„ H (^) — ^(I^PB -ffo(^)l^) along special directions related to the spec- 
tral decomposition of matrix Eq. Here our approach is different and based 
on wavelet technics. 

5 Complements and proofs 

5.1 Role of the real diagonalizable part of the anisotropy E 
of the analysing spaces Bp^^(M.'^,E) 

We will first prove that measuring smoothness in the general Besov spaces 
BpgiM.'^^E) may be deduced from the special case where the matrix E is 
diagonalizable. To this end, we show the following embedding property: 

Proposition 5.1 Assume that Ei G 8'^ and E2 € £^ have the same real 
diagonalizable part D. Let \ ■ \t^^ (resp \ ■ \t^^ )) a {W^, ^E^ ) (resp (W^, *E'2 ) ) 
pseudo-norm. Then for any a > and any {p,q) G [l,-|-oo]^ one has, 

, ^{R^Ei) ^ B'^^^{W',E2) ^ B" , ,,(M^^i) . 

p,q,[log| PminCO) ' p,lJ,|loglPmm{-D) 

(5.9) 

As a direct consequence, we obtain Corollarv 15.11 

Corollary 5.1 The anisotropic local critical exponent 

ax,iociE,P,q) = sup{s > 0, X{-) G B;^^i^^{R'',E)} , 

of any Gaussian field {X{x)}^^^d in anisotropic Besov spaces Bp^g{W^,E) 
depends only on the real diagonalizable part of E. 

Note that this result does not depend on the studied Gaussian field but of 
the analyzing functional spaces. Hence, it does not give any information 
about the anisotropic properties of the field. 

We now show Proposition 15.11 The proof of this result relies on the 
following lemma : 
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Lemma 5.1 Assume that Ei and E2 are two matrices of £^ having the 
same real diagonalizable part D. Then there exists two positive constants c\ 
and C2 such that, for all x G M*^, 



Ci\x\tE,^il + \logi\x\tEj\) ".ninC^) < \x\tE^ <C2|xh^^(l+|log(|x|ts^)|)''n.i„(^) . 

(5.10) 

Proof of Lemma 15.11 Using polar coordinates associated to ^E^, one has, 
for X G R'^, 

x = r^^@,{r,@)emlxSf' {\-\tE^). 
Denote Fi = Ei - D, F2 = E2 - D. Then 

< r\r~^^'^r^^Q\tp , 

because Fi,F2,D are pairwise commuting matrices. Observe now that -Fi, 
have only pure imaginary eigenvalues. Hence, by Lemma 2.1 of [.7J, one 
deduces that for any e > 

\x\tE^ < Crma.-x{\r~'^2r'^lQ\Pmin(0)-^ ^ |7.~*^27.'P'l9|pmax(0) + E:) , 

where | ■ | denotes the usual Euclidean norm. Denote || • || an operator norm 

on Mrf(R). Since Q belongs to the anisotropic sphere Sq ^ (| • |t_E^ ) which is 
compact, one has 

\x\tE < Cr ma.x{\\r~^^2r^^ \\ Pmin(D)-e ^ ||7--'-'^2 7-*-^l || Pmax{D) + e ) 



< Cr (1 + I log(r) I) "minC-D)-- 

< Cr (1 + |log(r)|)fmin(o) , 

for e > sufficiently small. We then proved Lemma 15.11 We now show 
Proposition 15. li 

Proof of Proposition 15. IL Using two anisotropic Littlewood-Paley anal- 
ysis associated respectively to matrices Ei,E2 and D and the lemma above, 
we deduce (15. 9p . Indeed, for any i G {1,2}, let (V'^')jgN an anisotropic 
Littlewood-Paley analysis of Besov spaces B^g{R'^,Ei). By definition, 

supplV?-) C {e G 1 < \C\tE^ < 4} , 

for i G {1,2}. Then there exists some jo G N such that for any j G N, one 

has 

supp(^f) C {C,2^-i<|ehs^ <2^+i} 



-' -"J o ■ (D) 
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Hence 



\ 



V 



Define the conjugate of q, that is the positive real satisfying 1/q+l/q' = 1. 
The last inequality and Cauchy-Schwartz inequality imply that for some 
C>Q, 



( 



ii/*V'fiiip < c-aog^jr^^' 



E 



J JO . /jj-v 



< C7(log2i)^/^'||Vo^^||ii 



^min'-^-' 



E 



LP 



(=j-jo- 



Then we can give the following upper bound of Yl 2^^"^ 11/ * V'f^ Wlp' 



i=i 



^'minv-^/ 



IK/'V'f)!!'' 



LP 



dlog2(0 



E (iog2jr/''2^''' 



V 



-JO- 



Pmin(^) 



T I I d\og2(J) 

< C E II/* Vf'llip2^'^^'^'''''"'"^-°^+^ . 



Let now J tends to oo. It yields the embedding 



B. 



p,?,|log| Pmin(O) 

Permuting E'l and yields the other inclusion. 



5.2 Local regularity in anisotropic Besov spaces of the stud- 
ied field 

In the previous section, we proved that we can restrict our study to di- 
agonal Besov spaces. This point is crucial for the proof of the regularity 
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results stated in Section HI Indeed it allows us to use tools that are only 
defined in the diagonal case, as anisotropic multi-resolution analysis and 
anisotropic wavelet bases. The aim of the following subsection is to recall 
the constructions of these wavelet bases. 



5.2.1 Orthonormal Wavelet bases of (diagonal) anisotropic spaces 

In this section, we assume that the anisotropy D of the analyzing space is 
diagonal (with positive eigenvalues), namely that 



D 



Ai 



In addition we also assume that Tr(L') = d. Our main tool will be anisotropic 
multi-resolution analyses defined by Triebel in [36j. 

Let {Vj,j > 0} be a one-dimensional multi-resolution analysis of L^(M). 
Denote by ^/^^ (resp. tp^^) the corresponding scaling function (resp. wavelet 
function). 



Notation 5.1 We denote by {F,MY* the set 

{F,mY* ={F,mY\{{F,--- ,F)} . 
For j G N, we define the set P {D) of {F, AI}'^ x N'^ in the following way. 
. Ifi = 0, /0(Z)) = {((F,... ,F),(0,--- ,0))}. 



If j > 1, P{D) is the set of ah the elements (G, 7) with G G {F, M^* 
and 7 G N"^ such that for any r G {1, • • • ,d} : 

If Gr = F,^r = iU -l)Xr], 

If Gr = M, [{j - l)Xr] < 7r < [j K] • 

Finally, for j G N and {G,j) G P{D), we will denote by -Dj,G,7 the 
matrix defined by 



D 







Id) 



Finally, let us define the family of wavelets as follows. For j G N, (G, 7) G 
P{p) and k G Z"^, we set 

*^,G,,(^) = (V'(''^)(2''-«-3;-fe), 
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with 

jp(G) = . . . . 

The anisotropic wavelet bases yield a wavelet characterisation of anisotropic 
Besov spaces (see |33 and [36], Theorem 5.23). 

Theorem 5.1 

1. The family 1^2^^^^^^^^^^^^, j G N, {G,j) G P{D), k G Z^'j is an or- 
thonormal basis of L'^iW^). 

2. Let {^k^''^)j^n,(G,'y)£p{D),keZ'' the family constructed from tpp and 
ipM Daubechies wavelets with, for some u G N, 

Let < p,q < oo and s,A^ G M. There exists an integer u{s,p,D) 
such that if u > u{s,p,D), for any tempered distribution f the two 
following assertions are equivalent 

(b) / = ES^,G,7^k7 

i,G,7 V k 

the convergence being in S'(M.'^). 
The above expansion is then unique and 

cla,,=<f,2^'^''-''-^^lG,,> ■ (5.11) 

Remark 5.1 An analogous result is stated (see [36], Theorem 5.24) replac- 
ing Daubechies wavelets by Meyer wavelets. In that case, u = +oo. 

We now prove our regularity results about the sample path of {Xp^^^Ho{x)}xeR'i 
based on wavelet characterization of Besov spaces. 

5.2.2 Local regularity of the field {XEo,Ho{x)}^^^d in anisotropic 
Besov spaces B^^g{W^,Do) 

Assume that we are given a Gaussian field {XEo,Ho{^)}xeR'i of form ()2.3p 
where Eq G SJ^ and Hq G (0, pmin(^'o))- The aim of this section is to prove 



< +00 
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Proposition 5.2 Let 1 <p,q< +oo. Define 5 on (0, +oo] as follows : 

6(p) = I ifP = +oo, 
\ 1 otherwise. 

Then one has 

1. For any j3 > 1/q + (i//9min(-E'o) + ^{p)j almost surely, the sample path 
o/{^PBo.^o(^)}xGR<' belongs to 3^°^^^^^^^^ ^^^JW^ , Dq) , 

2. For j3 = 1/q + d//Omin(-£'o) + ^{p), almost surely, the sample path of 

p,q,|log| ^ ,loc^ 



{XpEo,Hoix)}xm'' does not belong to , ^^^^(R'^,L»o) 



Adapting to our setting a result of [22], we first remark that there exists 
C°°(M'^\{0}) (lR'^,^o) pseudo-norms 

Lemma 5.2 Let Eq S and ip he a C°° function compactly supported in 
\ {0}. The function p defined on , by 



p{x) = I (p{a °x)da , 



is a (M , -Eq) pseudo-norm belonging to C^{W \ {0}). 

We now prove that the belongness of the sample paths to anisotropic Besov 
spaces of any OSRGF of the form {-'^^p./ZolaeRd do not depend on the (M*^, ^Eq ) 
pseudo-norm p involved in the construction of the field. 

Lemma 5.3 Let Eq € and pi,p2 two ( £^q,M'') pseudo-norms. Denote 
respectively {Xi{x)}^^^d and {Xi{x)}r^^^d the two OSSRGF defined from 
pi and p2. Then, for any s > 0, /? E M, {p,q) G (0,oo]^, a.s. Xi belongs to 
^P,g,|iog|/5(^'''^o) iff a.s. X2 belongs to Dq)- 



Proof. Remark first that using the same approach than in Lemma 2 of j8] 
and an anisotropic version of Kolmogorov Centsov Theorem we can prove 
that a.s. 



X I-)- 



j (ei<-''«>-l)p-^'>-'^/2(e)diy(e),x^ / (e^<-'«>-l)p-''"~'/'(OdH^(e) 



both belong to C^{K,Do) Bpg{R'^,Do) for any compact subset K of 



M and any r > min;j^g5p(^) A. Since any (^EqjR ) pseudo-norms pi, P2 are 
equivalent, Lemma [5. 3 1 is then a straightforward consequence of Theorem 1.1 
of |12) applied with B = Bp^{W^,Do) which is either a separable Banach 
space either the dual of the separable space B = Bp,^^,{W^, Dq) with p',q' 
the respective conjugates of p, q and to 

fx = Pl'''"-\\>R and fy = P2''"-%\>R , 
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and 

fx = pf'''-\\>R and fy = Pi'''-\i>R , 

successively. 

Thus, using Lemmas 15.21 and 15.31 we assume without loss of general- 
ity from now that the (M'^,£'o) pseudo-norm peq, used to define the field 
{Xp^^^Ho{x)}x£m'i belongs to C°°(R'^ \ {0}). We shall use this assumption 
when proving that the wavelet coefficients of the field are weakly dependent 
(see Section [6]). 

Observe that, to prove our local regularity results, we have to investigate 
the sample paths properties of ipX for any function ip G that is for 

any /cq G Z*^, tq > and any ip G 'D{M.'^) satisfying supp((/?) C B£)^^{kQ,rQ). 
Since the Besov spaces are invariant by translations and dilatations, we 
may assume that ko = and ro = 1. We have then to study the sample 
paths properties of the field (pX for any function ip G P(M'^) such that 
supp(v7) C -Bdo(0, 1). 

Our results come from the series expansion of Xp^^^Ho in a Daubechies 
anisotropic wavelet basis (see Section [5.2.11 iust above). Recall that for any 
j G N, (G, 7) G Ij{D), the wavelet coefficients of Xp^^^H^ are defined as 

Fix now I • a (M°',Z)o) pseudo-norm. Define Tq{Dq) = for j = and 
for any j > 1 

r,(Z)o) = {fc G \k\D,<i2^]. (5.12) 



Thereafter set 



and 



We will investigate separately the local sample path properties in anisotropic 
Besov spaces of the two Gaussian fields X^^^ and X^^^ . We first prove 
that 

Proposition 5.3 Let I < p, q < 00. 

1. Almost surely, for any (3 > 1/q + (i/pmm(-E'o) + ^(p), the sample path 
of the field {Xl^l^j,^{x)},^^, belongs to B^;^^^^^^^,{R^, Do). 
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2. Let if such that supp((^) C Bdq{0, 1) and satisfying 

^ = 1 on (0,1/2) . 

Then almost surely, for /? = 1/q + d/puimiEo) + ^{p) the sample path 
of the field {^jX^^^^ j^y(x)}^gRd does not belong to B^'^^^^^^^_f^{R'^,Do). 

Proof. The proof uses several technics introduced in [11]. The result comes 



from a comparison between X^fcer, 



{Do) IS-,G,7l 



1/p 



and 



from Lemma l6 . 1 1 which gives an estimate of 



1/2 



Set 



1/2 



and 



1/2 



(5.15) 



for any j G N, (G,7) G Ij and A; S rj(L'o)- We need to distinguish two 
cases : p ^ oo and p = oo. In each case, we prove successively points (i) 
and (ii). 

Assume first that p ^ oo and let us prove point (i) in this case. The 
definition of the sequence (gjc^) and the stationarity for any {j,G,j) of 
the sequence {Cj q G Z'^) implies that for any j, G, 7 



i/p 



E 



2m1/2 




Use now the weak correlation of the wavelet coefficients and the two esti- 

1/2 

and of Uj = card(rj(Z)o)) respectively proved in 



mates of 



IE(|c?^j2) 



Lemmas 16.11 and 16.21 One deduces that the following inequality holds for 
any j > : 



yk&j(Do) 



< C2^(|-^«)j'^* 



\9j,Ga\ 



k£Vj{Do) 



(5.16) 



where 



d d 
d* = I 

2pmin(-E'o) P 



Lemma 16.31 stating a central limit theorem for the sequence (gjQ^) and 
inequality (j5.16p then prove point (i) of the proposition for the case p < 00. 
We now prove point (ii) for p ^ 00. Set now 



r' (Do) = {A: G < 2Vj} . 
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Using the assumptions on the support of remark that for j sufficiently 
large and for any k G T'j{Do), one has 

Use the same arguments as in the proof of Lemma 16.21 and deduce that 
n'j = Card(r^(Z)o)) ~ Since T'^{Dq) C rj(L»o), a similar approach 

to above then yields that for some C > and for any j > 1 




which directly implies point (ii) of the proposition. 

If p = oo, a similar approach implies that almost surely there exists some 
Ci, C2 > such that for any j, G, 7 



Ci2-^-^0j-'^* [ J—- sup \g^^aM < ( E IS'C^I" 



and 



E |cf,^,,r < C22-^'''^r sup \gia,,\ 



with 



d* 



2pmin(-E'o) 

Lemma 16.41 and the inequality just above then implies the result stated 
in point (i) for the case p = 00. The proof of point (ii) for p = 00 also 
follows from the above inequality replacing Tj{Do) with r^(L'o) as in the 
case p ^ 00. 

(2) 

We now investigate the sample paths properties of {X^J Ho(^)}a;GR'*- 

(2) 

Proposition 5.4 Almost surely, the sample path of the field {X^^ Ho{x)}xgb 



belong to B^'^ i^^{M.'^,Eq) for any 

0<Ho<H' < Pmm(-Do) = /9min(-Eo) , 

and any 1 < p,q < 00. 

Proof. Using the transference results of [36] (see Theorem 5.28) and the 
usual embedding of isotropic Besov spaces defined on bounded domains one 
remarks that 

Cf+^(M^I)o)ci?^,,,,„,(M^Z?o), 
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for any 1 < p,q < oo and any s,e > 0. It then suffices to prove the result 
for p = q = oo. 

Let now consider H' € (i^o, Pmin(-E'o)), e > and (f G T>{M.'^). Recah 
that we assumed that 

supp((^) C BooiO, 1) = {x, \x\do < 1} , 

and < 99 < 1 on M'^. We denote by Y the random field ^X^p^ 

We will give an upper bound of \Y{x + h) — Y{x)\ for any given x in 
i?£)p(0, 1) and h sufficiently small. Observe that 

Y{x + h)- Y{x) = Yi{x, h) + Y2{x, h) , 

with 



^7W 



Yiix, h)=Y, Yl 4g,M^ + ^) - V'(^))^i,G,- 

j,G,jk^rj{Do) 

We first bound Yi{x,h). Let e = 1 — H' / pminiEo) . We now use that 
If € ^^^ioc(^'^)' ^ sufficiently small and x belongs to the compact set 
B£)q{0, 1). Hence , by Lemma [6. 41 and the fast decay of the wavelets, almost 
surely for any M > and for some C > one has 



1 



+ \k- 2^^^G,-rx\)M 

Here we denoted d* = 1/2 + (i/pmin(-Eo). Further, by assumption on k and 

IA:|a, > i2^' > i|2^^'«'-x|Do . 

Since | • \do is a (M'^,Do) pseudo-norm, by the triangular inequality (|3.5p . 
one deduces that for j sufficiently large 

\k - 2^^'^'^x\do > C\k\Do ■ 

for some C G (0,1). Then by comparison between j • {jj^ and the usual 
Eucidean norm, one deduces that there exists some a > such that for j 
sufficiently large and any x in Bdo{0, 1) 



\k - 2^J'G,7x| > {\k\/2Y 
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Then 



Since, for M sufficiently large 

(l + l^l") 



i,G,7 fc^rj(Do) 



one has almost surely |Y'i(x,/i)| < C'\h\j!)^^. 

By the same approach, we can bound Y2{x,h). Indeed, using the fact 
that 99 is bounded and the mean value theorem for "^jc-y^ then prove 
that almost surely for some C > 

|y2(x,/i)i < c ^ /*2-^-^«i2^^.G,7/i| sup r 



+ \k- 2^i'G,jy\)M 



The end of the proof is exactly the same as above remarking that 

for some 5 > 0. Proposition 15.21 then follows directly from Propositions 
and EH 



5.3 Proof of regularity results in anisotropic Besov spaces 
with an anisotropy commuting with this of the field 

The following proposition extends the results of Proposition l5.2l in anisotropic 
Besov spaces Bp^^{W^, E) with E £ commuting with Eq. 

Proposition 5.5 Let 1 < p,q < +00, e > and E £ £^ commuting with 
Eq . Then 

1. Almost surely the sample path of {Xp^^^Ho{x)}x<m.''^ belongs to ^ 

2. Almost surely the sample path of {Xp^^^Ho{x)}.^^^d does not belong to 

The proof is made in several steps. First we need to compare Besov spaces 
with different commuting anisotropies. 
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5.3.1 A comparison result between Besov spaces with different 
commuting anisotropics 

Since E and Eq are commuting, we can then assume (up to a change of 
basis) that Dq and D are two diagonal matrices of the form : 



with 



\ 




D 



\ 







A5k 



A? 



(5.17) 



(5.18) 



Proposition 5.6 The notations and assumptions are as above. For any 
a > 0, /? € M and p,q (0, +oo], one has the following embedding 



p,q,\ log ^ 



B 



P,<?,|log 



The proof is straightforward and based on finite differences characterization 
of Besov spaces given in Theorem 5.8 (ii) of [36]. 

5.3.2 Proof of Proposition [5T5l 

We only prove the second point of Proposition 15.51 since the first one is a 
straigthforward consequence of Propositions 15.11 and 15. 6[ To this end we use 
the following characterization of anisotropic Besov spaces B%W^,/^) with 
diagonal anisotropy A (see Theorem 5.8 of [36]) : 

Proposition 5.7 Let A a matrix belonging to of the form 



A 



(ax 
a2 

V 



0\ 


Oid) 



and s G (0,pmin(A)), = [s/ae\ + 1 for any i = I 
B^^p{R'^, A) if and only if 

d 

\LP 



, d. Then f G 



where (e^) is the canonical basis o_ 
iAlf)ix) = fix + h)-f{x 



and where as usual, for any x, /i G 
, (Af /)(x) = (AiA*^-V)(x) . 
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Remark 5.2 If for ^ = 1, • • • , d,s G (0, ae) then Me = 1. 

Proof. This statement is proved in Theorem 5.8 in |36j . 

We now prove Proposition 15.51 We first remark that we have only to 
consider the case where Eq = Dq. Indeed, let | • l^;^ (resp | • \do) be a 
(K'^,Eq) (resp a (M'^,L'o)) pseudo-norm. Lemma l5.ll then implies that for 
any e > and any |^| sufficiently large 

Hence Theorem 1.1 of fl^ applied successively with fx = l|f|<ij|'?l£o^°~'^i Jy = 

above and Proposition 15.51 proved in the case Eq = Dq yield the result in 
the general case Eq £ . 

From now, we then assume that £"0 = -Do and that the (M*^, -Dq) pseudo- 
norm involved in the construction of the studied field belongs to C°°(M'^\{0}) 
which ensures the weak correlation of the wavelet coefficients. 

As in the proof of Proposition 15.21 we use an expansion of the Gaussian 
field Xp^^^Ho in a Daubechies wavelet basis and we define, 

j,G,7fcerj(Do) 

and 

X^^\^)=Y. E 4G,7('^)^k7(^) • (5-20) 
i.G:7fe0rj(Do) 

where fj{DQ) = {A; G Z'^, |A;|do < Ci2^. 

As in the proof of Proposition 15.21 we see that, for Ci sufficiently large, 
almost surely X^^) belongs to ^^^^^'/^^'^^(K'^, D) for any 1 < p, g < cx) and 
e > sufficiently small. 

We then have to prove our a.s. non local regularity results for the Gaus- 
sian field Remark now that since the multiresolution analysis is com- 
pactly supported so is X'^^\ To show point (ii) of Proposition 15.51 it is then 
sufficient to prove that a.s. the sample paths of X*-^-* does not belong to 
Q\r^Ho/\t+e{^d^ D) for any e > 0. 

Set M = [s/\rri\ + 1 which may be greater than one. In view of Propo- 
sition 15.71 we shall then give an almost sure lower bound of 

ll(A,t^(^))(x)||i,= / |(A,^f^X«)(x)Rx, 

for any p > 1 and any t of the form t = 2~[-^"^'"] where jo is a fixed non- 
negative integer. 
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Set 



/boA?] 



A 



JO 



V [joA^ 

Observe that if t = 2~^^°"^'^^, one has te^ = 2~^ned- Remark also that 
X^^^ can be written as the sum of its low frequency component and its high 
frequency component, namely that 



with 



» and X'^l,{x) 



E E E 



i<io G,7 fcgr^.(£)g) 
Using the triangular inequality, one has 

||(Aff^xW)(x)||L. > \\{Afi^x'^l){x)\\L. - \\{Af^^X^ll){x)\\L. (5.21) 
To give a lower bound of ||(Ai^^^X(i))( 

IIlp ' shall then give a lower bound 
of \\{A^lX^j^l){x)\\% and an upper bound of ^(Ai^^xS)(x)||i,. 

Let us first give an upper bound of ||(A4^^X^^^)(x)||£^p. We suppose that 
the multiresolution analysis is s smooth for some s G (i^o, /Omin(-Do)). By 
the finite differences definition of the spaces B^^^{W^, Do) and the fact that 
for any M > 1 |Af < J^eU \f{x + ih)-f{x + i£-l)h)\, one has 



EE E 4c.(A«.„^^*i,„)M 

j<jo G,7 fcgrj(Do) 



LP 



j<jo 



E E s-g,7^'g. 

G,7 fcer,(A)) 



i?|,^(IR<*,Do) 

Use now the wavelet characterization of the homogeneous Besov spaces 
B^^oci'^'^, Do)- Then for some C > 

i/p 



E E S-,G,7^j,G,7 
G,7fcerj(Do) 



<c2^-(^-^/rtK^ E ls-G.71 



,G,7 fcerj(Do) 



As the proof of Proposition 15.21 we can estimate a.s. "^f^^f .^j^^-^ \cjq^\'p. 
Hence, we deduce that there exists an a.s. positive constant C such that 



LP 



j<jo 



< 02'^°" ^ 2^'^^*"^^ • ( 2^Pj^Pmin(EQ)2-i^0 ] < (J' 



(5.22) 
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where C is not depending on jo nor s. 

We now give a lower bound of ||(Ai^^^X^^j,)( To this end perform 

the change of variable x = 2~^ny and deduce that 

ll(A,t^S!.)(^)irL. = 2-Tr(A,,)||(^M )(2-^^oy)f^^ (5.23) 

By definition of X^^p one has 

i>jo+l G,7 fcgrj(Do) 

Define now for any j' > 1, any G G {F, Af j'^' and any 7' € N*^ such that 

[{f - l)Xr] - 2 < 7; < [j'K] + 2 
the family of functions 



where 



/7i 

■•• 
Vo ... 



7™/ 



and for / = and any k ^ 11^ h^-. 



•i',(F,..,F),(0,...,0)(^) = 

Observe that this is a family of inhomogeneous smooth analysis molecules 
in the sense of Definition 5.3 in |9j. Further, if j' = j — jo and 7 = ([jo^^rDr- 
one has 



x(;j,(2^^^o (y + e,)) - (2-^^oy) = E E E ^',G,y^' ,G,y (y) , 

/>1G,7' fcgrj(Do) 

with = 4+^.^^^ .^,+- if (G,y + 7) € P'+^°(Do) and ^,7 = « ot^er- 

wise. Hence 



(1) 



2 TO Ed 



2 Wed 



IfiO ,, ,(K<*,Do) 



> 



We use once more an a.s. estimate of ^^i. |cj(^^|p as in the proof of Propo- 
sition [521 Since as in Lemma [6.21 we can prove that Card(rj/+jo(Do) ^ 
2O +jo)d^ Hence there exists an a.s. positive constant C such that 



E 



Ufc IP > ')-(]' +3o){Hop-d) 

IS''+io,G,7l - ^ 
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Hence 

2 ■'Oed 



Use now the last inequality and relation (I5.23p . Then a.s. 



We deduce that a.s 

-1 

+00 „2-bo^ml 







io=o 

+ 00 



"l"°° ■ -sp/X„ 



> C7^ j-<i2--'«-f^«p(2-b"^-l) ^2-[^«^- 
io=o 

+ 00 

io=o 

Since 

io=o 

if sX^/Xm — Hqp > 0, we deduce that a.s. 
Jo 

for s > HoXm/X^. Using (|5.22|) and the triangular inequality (|5.21|) . it ends 
the proof of Proposition 15.51 It also implies directly Theorem 14.11 

6 Technical lemmas 

Our results about smoothness of the sample path are based on the following 
lemma 

Lemma 6.1 Assume that the anisotropic multi-resolution analysis consid- 
ered is and admits at least one vanishing moment. 

Let {Xp^^^Hg{x)}rc^^d the Gaussian field defined by 112. 3\) with p = peq. 
Assume also that the pseudo-norm involved in the construction of this 
field is at least C^(]R^ \ {0}). Then the wavelet coefficients of the random 
field {Xp^^^Ho{x)}x£Rd are weakly dependent in the following sense 
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1. There exists some Co > such for any j > 1, (G, 7) € Ij and {k,k') E 

|E(c,^,G,,S^,G,.)l < go \^|fe_^.j • (6.24) 

There exists some Ci, C2 > suc/i i/iai /or any j > 1, {G,j) € /j and 
any A; G (Z'^) 

C7^j-rf/Pmin(So)2-2i/^o < E(|c^^G,7p) - C2i'^/^"""(-^°)2-2j'^o , (6.25) 

Remark 6.1 Theorem 1.1 of [?] imply that, studying the sample paths 
properties of OSSRGF, we can always assume that the pseudo-norm peq in- 
volved in the construction of this field is C°°(]R'^\{0}). Then the assumptions 
of Lemma l6. II are satisfied. 

Proof.Since the anisotropic multiresolution analysis admits at least one 
vanishing moment, one has V'm(O) = 0. Further, for any j > 1, (G, 7) G Ij 
and for all k € Z'^ 

This formula implies that (set C = 2~ ^) 

bmce using Lemma 15.11 and the inequalities 

Ci2^ < |2'^^.G.7 (1^,, < C22^ imply that 

E(|c^,G,7l') ^ Ci2-2^-(^o+'^)2^-'^ / |V^(^)(C)PlCl7^'''°"" (i+log(|C|n,)+Jr'/''"""^''°^ dC , 

and 

E(|c^,G,7l') < C22-2^-(^«+'^)2^-'^ / IV^^^HOI'ICC""" (l+log([Cb„ )+jf /^-'"(^«)dC . 
We then proved inequalities (|6.25|) . 

To prove inequality (|6.24p remark that for any G {1, • • • , d} 

Set C = 2 ^''^'T C and integrate by parts with respect to C^. Hence 
{h-k'^MclcAGn) = 2^'^'-(''-«-) / {h-k[)e<''-^')^\^^^\0\^pE,{2''^^o.-, C)-2^^o-rfdC. 
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Recall that the pseudo-norm may be assumed to be C°°(R'^ \ {0}). Since 
the multi resolution analysis is 



An approach similar to the proof of inequalities (j6.25p yields 



2iTr{D,,G,^) 



dC 



(6.26) 



Further, differentiate the homogeneity relationship satisfied by pt^^ and 
deduce that for any a > and z 



a^o (grad(pt^^ ))(a^o z) = a(grad(pts^ ))(z) 



(6.27) 



For any y € M'^ \ {0}, let r = \y\t^^ ■ Then set j = [log2(r)] and remark that 
= \2~^*^o y\t^ £ [1/2,2] and hence that belongs to the compact 
set C(l/2,2, %) = {e, \e\tE^ G [1/2,2]}. Relationship (fOT]) applied with 
a = 2^ and z = @ then implies 



2^^.o,. grad(p5,^))(2^*^o G) = 2-^"^o + 45..G,.2i(grad(/>5,J (9)) 

Take the norm of each member of the equality and deduce that for any 
y eM.'^\ {0} satisfying j = [loga |y| tijj 

|2^^-.G,.^(p,^ ))(y)| < C72J-|2-^"^o+5?,,G,, I 



where C = supQ,^c{i/2,2,%) k^^^^ip^) (0)1- 
Lemma 2.1 of [6] and the definition of j imply that 



|2^..G.. gvad{ptE^)){y)\ < C2^[jf/''- < C\y\tE^ |log(|yh^^ )|'^/Pn.in(i?o) . 
Set now y = 2 ^. One has 

|2*^.,G.. {^d{p^^ )){2^..o.-, Q\ < CI2*^.,G,. Cl*^^ |log(|2'^.,G,. Cl*^^ )|'^M.in(i^o) 

< C2^-(j+|log(|Ch^J)|)2'^/''-"(^«)|C|5;^. 

Since for any ^ G {1, • • • , d} 



27< 



'j,G,7 C) 



< 



2^..G.7(grad(/9ijJ)(2".>G,. ^) 
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it yields the following inequality 

-m2 / n \ 

< 2-2jHo]^j^2d/p^in{Eo) 



ei(^-fe')C|^(G)(C)|2 ^^^^ Q 

(6.28) 

Combining inequalities (|6.26p and ()6.28p then yield inequality ()6.24p . 
Remark now that 

Lemma 6.2 Let Dq an admissible diagonal anisotropy satisfying Tr(L'o) = 
d. Recall that Tj^Dq) is defined by ^5.1^) . There exists some Ci,C2 > 
such that 

Ci/2-''^<card(r,(Z)o))</2^''. 

Proof. Indeed, since the norms | • and | • |^oo on M*^ are equivalent, there 
exists some Ci , C2 > such that 

Ci max Ikf^/^' < I^IA) < Camaxlfc^l^/^* . 
The conclusion follows since it is quite clear since that 

card{fc,max|/c£|i/^^ < j2^} = card{A;, max |/c£| < = JJ (^j^^2^'^^) = /2^ , 

using the fact that Xi + ■ ■ ■ + = d. 

The proof of Proposition 15.31 is then based on the two following results 
which are a slight modification of Theorem ILl and IL7 of We recall 
the proofs for completeness. 
We denote 

cp = n\9lGa\')- 

where {g^Q^) is the stationary Gaussian sequence of the normalized wavelet 
coefficients defined by (I5.15p . Since, by Lemma [6T| the wavelet coefficients 
are weakly dependent, we can state a central limit theorem for the sequence 
(ffjG 7)ieN {G'y)&i- feer which is a slight modified version of Lemma IL4 of 



Lemma 6.3 Let p G (l,+oo) and {g^Q^) the Gaussian sequence defined 
by ^5.15\) . Set Uj = CaidiV j{DQ)) . Then almost surely when j — > 00 



fcerj(Do) 
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Proof. By Lemma fc.ll the sequence (gjQ^) is weakly correlated in the sense 
of |ll]-that is satisfies the assumption (H) of pT]. We follow the main line 
of [H] and first give an upper bound of 



E 



Using the same approach that in [TT] (see Lemma II. 3) we get that 



E 



(fc,A:')Gr2 



with Cj = j2'^/Pmin(£^o) i^y -vveak correlation of the wavelet coefficients. Set 
i = k — k' . Hence 



^ (l + \k-k'\)^-^ E (1 + 1^1)2 - ^^'-^'^^ E (1 + 



\2 ■ 



Remark now that 

with 6 = 2/pmax(-E'o) > by comparison with an integral and Proposi- 
tion 2.3 in [7J. 

Thereafter the end of the proof is exactly the same that in Theorem II. 1 
in [n]- 



In an analogous way, one can give a result on the asymptotic behavior 



of 



^=^= (max I G -7 1 I 

.g(n,)| V^er,'^'' '^7 



Lemma 6.4 Almost surely 
< liminf 

j->oo 



max \gj^G,j\ < limsup 
n^)\ VfeGr, J j^oo 



max\gj^G,-y\ < oo 



Proof. The proof is exactly the same than these of Lemmas II. 8 and 11.10 
in 
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